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SYMMETRIC  FRACTALS  GENERATED  BY  CELLULAR 

AUTOMATA 
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Kardinaal  Mercierlaan  94,  3001  Heverlee,  Belgium 
E-mail:  {Andre. Barbejvanhaes} Qesat.kuleuven.ac.be 

We  demonstrate  how  to  construct  fractals  which  are  generated  by  a combination 
of  a cellular  automaton  and  a substitution.  Moreover,  if  the  substitution  and  the 
cellular  automaton  exhibit  certain  symmetry  features,  the  fractal  will  inherit  these 
symmetries. 
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1 Introduction 

It  is  a well  known  fact  that  certain  cellular  automata  can  generate  fractals7, 8,9,14,15 . 
These  fractals  can  be  described  in  terms  of  hierarchical  iterated  function  systems12 , 
graph  directed  constructions11  or  as  a mixed  self-similar  set2. 

In  this  note  we  consider  cellular  automata  which  have  certain  symmetries.  After 
the  introduction  of  some  basic  concepts  and  formalism  about  substitutions,  fractals 
and  cellular  automata  in  Sections  2 and  3,  we  shall  show  in  Section  4 how  to 
construct  fractals  generated  by  a cellular  automaton  such  that  the  fractal  inherits 
the  symmetries  of  the  cellular  automaton. 

2 Substitutions  and  fractals 

In  this  section  we  introduce  the  concept  of  a substitution  which  is  needed  for  our  ap- 
plications. We  discuss  how  substitutions  generate  a compact  subset  which  exhibits 
a hierarchical  self-similar  structure,  and  thus  is  usually  a fractal  set. 

We  will  consider  sequences  (fj)jez,  where  fj  belongs  to  a finite  commutative 
ring  1Z  with  0 and  1.  It  will  be  useful  to  represent  this  sequence  as  a formal  Laurent 
series  with  coefficients  in  7 Z,  as  / or  as 

l = f(X)  = Y/fjXj- 

je  z 

The  set  of  all  such  Laurent  series  will  be  denoted  by  1Z{X).  The  support  supp(/) 
of  / £ 'R-(X)  is  the  set  {j  \ fj  0}.  / is  called  a Laurent  polynomial  if  supp(/) 
is  finite.  The  set  of  all  Laurent  polynomials  is  denoted  by  1ZC(X).  If  supp(/)  C N 
and  is  finite,  then  / is  called  a polynomial.  The  set  of  polynomials  with  coefficients 
in  TZ  is  denoted  by  1Z[X). 

The  set  of  all  maps  a-.IZ^IZ  equipped  with  the  addition  and  composition  is 
denoted  by  Abb(7Z).  With  0 we  denote  the  map  r 4 0 and  with  1 we  denote  the 
identity  map  r4r.  The  subset  of  all  maps  a : TZ  -4  TZ  with  a(0)  = 0 is  denoted 
as  Abb0(7£).  If  a G Abb  {TZ),  then  a induces  a map,  also  denoted  by  a,  from  7 Z(X) 
to  1Z(X)  which  is  defined  as  a(f)(X)  = ■ 

For  k € N the  set  {0, . . . , k — 1}  is  denoted  by  [k]. 
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A ^-substitution  transforms  a sequence  (fj)je z into  a new  sequence  by  replacing 
each  element  fj  by  a string  of  k elements  £o(/j)£i(/j)  ■ ■ -fk-i(fj),  where  £,  6 
Abb  (72).  Here  follows  the  formal 

Definition  2.1  A fc-substitution  £ is  a k-tuple  £ = (£;)ie[fc]  £ (Abb(72))*  which 
defines  a map  £ on  TZ(X ) given  &«/ 

fc-i 

£(/)  = £ *'&(/)(**), 

J=0 

where  fj(f)(X)  = &(/.?)  AP  #£  £ (Abbo(£>))*,  he.,  when  0 is  replaced  by  a 

string  of  k 0 ’s,  then  £ is  called  a regular  k-substitution. 

Remark  1.  A /’-substitution  £ = (£0;e[fc]  is  nlso  defined  by  a polynomial  T\{X)  — 
g Abb (72.) [A']  of  degree  less  than  k and  we  write  £(/)  = P^(X)  o f(Xk), 
where  the  product  is  the  product  induced  by  the  composition  of  maps  and  the 
coefficients  fj  of  f(X)  are  considered  as  constant  maps  rn/,.  The  polynomial 
P^  is  called  the  substitution  polynomial 

2. If  Q (X)  £ 72[X]  is  a polynomial  of  degree  less  than  or  equal  to  k — 1,  then 
£q(/)  = Q(X)f(Xk)  defines  a /c-substitution  with  £j(r)  = qjr. 

3. If  £i  and  £2  are  /c-substitutions,  respectively,  and  if  P5l  (X)  and  P^2  (X)  denote 
the  respective  substitution  polynomials,  then  £ = £i  o £2  is  a k2  substitution  with 
substitution  polynomial  P^(X)  — P^(X)  o P^2(Xk). 

We  also  need  to  consider  two-dimensional  sequences  (gi,j)i,je z with  Qij  £ 72. 
The  corresponding  formal  Laurent  series  representation  is 

g = f(X,Y)=  Y,gijXiYti 
i,jez 

TZ(X , Y)  is  the  set  of  all  such  two-dimensional  Laurent  series  (sequences).  Like  for 
the  one-dimensional  case,  we  say  that  g(X,  Y)  is  a Laurent  polynomial  if  gij  = 0 
almost  everywhere.  Like  for  the  one-dimensional  case,  two-dimensional  sequences 
can  be  transformed  by  a two-dimensional  ( k x /^-substitution,  which  replaces  each 
symbol  gitj  by  a (k  x Zc)-array  of  symbols  in  72.  Formally: 

Definition  2.2  Let  E = ( ] £ (Abb(72))fcxfc.  The  map  E : TZ(X,Y)  ->• 
1Z(X,Y)  defined  as 

k- 1 

E(g)=  £ XlYm^m(g)(Xk,Yk) 

l,m= 0 

is  called  a (k  x /c)-substitution.  It  is  called  regular  if  E £ (Abbo(72))fcxfc. 

Let  us  agree  here,  that  in  two-dimensional  sequences  and  arrays,  and  in  the 
corresponding  graphical  representations,  the  first  index,  which  is  possibly  associated 
to  the  symbol  X in  formal  Laurent  series,  increases  along  a horizontal  axis  which 
is  oriented  from  the  left  to  the  right;  while  the  second  index,  possibly  associated  to 
the  symbol  Y,  increases  along  a vertical  axis  which  is  oriented  from  top  to  bottom. 

Remark  1 . If  E is  a regular  ( k x /^-substitution,  then  H(72C(AT))  C 7lc(X).  In 
other  words,  the  set  of  Laurent  polynomials  is  invariant  under  regular  k-  or  ( k x k)- 
substitutions. 
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2. If  Hi  and  H2  are  (k  x k) -substitutions,  respectively,  and  PEl  , Pe2  the  respective 
substitution  polynomials,  then  Hi  oH2  is  a (k2  x A;2  ^substitution  with  substitution 
polynomial  Pei0e2(X,Y)  = PEl(X,Y)  o PE2(Xk ,Yk). 

In  order  to  draw  the  connection  between  fractals  and  substitutions  we  define  a 
graphical  representation  of  elements  in  7ZC(X)  or  in  TZC(X,Y).  We  only  need  to 
introduce  the  graphical  representation  for  elements  in  7ZC(X,  Y),  since  any  Laurent 
series  in  7L(X)  has  a natural  counterpart  in  1Z(X,Y). 

Let  (P(M2),d^r)  denote  the  set  of  non-empty  compact  subsets  of  R2 , where  dn 
is  the  Hausdorff  distance  induced  by  the  Euclidian  metric  (or  any  other  equivalent 
metric)  on  K2.  Since  M2  is  a complete  metric  space,  so  is  ('H,  dfj)5’6 ■ 

Definition  2.3  Let  I(i,j)  = [i,i  + 1]  x [j,j  + 1]  . Then  the  map  G : TZC{X,Y)  — » 
P(R2)  U {0}  defined  as 


G{g)  = { I(i,j ) 1 gij  # 0} 

is  called  graphical  representation. 

I(i,j ) in  the  definition  above  can  be  considered  as  a “highlighted  pixel”  at 
location  (i,j)  on  a graphical  display  which  coincides  with  the  E2 -plane  (as  agreed: 
y-axis  positively  oriented  downwards). 

Remark  l.G(g)  = 0 if  and  only  if  g = 0. 

2.  Any  map  7 : TZ  \ {0}  -»  P(R2)  induces  another  graphical  representation  G7  by 
defining  G 7(g)  = U{(i,  j)  + j(gi,j)  | 9i,j  # 0}.  This  corresponds  to  “ non-square 
pixels”,  where  the  pixel-shape  may  depend  on  gij. 

3.  The  concept  of  graphical  representation  is  quite  general:  when  S is  a ring  with  0 
and  1,  then  so  is  TZ  = SD,  with  Or  = (0, 0, 0, ... , 0)  and  1r  = (1, 1, 1, . . . , 1)  (both 
On  and  1-r  have  D components).  So  it  is  possible  to  speak  about  the  graphical 
representation  of  a two-dimensional  sequence  with  values  in  TZ  = SD , where  gij 
is  just  an  element  of  SD. 

The  connection  between  fractals  and  substitutions  is  provided  by  the  following 
theorem. 

Theorem  2.4  8 Let  H : 1Z(X,Y)  -»  1Z(X,Y)  be  a regular  ( k x k)-substitution  such 
that  for  all  r / 0,  there  are  i,j  such  that  fij(r)  ^ 0.  If  g 6 7 ZC(X,Y)  and  g^0, 
then  the  sequence 


n€N 


is  a Cauchy  sequence  in  ( TL,dH )• 

The  proof  relies  on  two  facts.  Firstly,  H being  regular  maps  Laurent  polynomials  on 
Laurent  polynomials.  Thus  the  above  sequence  is  indeed  in  7f(K2).  Secondly,  again 
due  to  the  regularity  of  H one  has  |G(H(^))  C G(.g)  for  all  Laurent  polynomials 
different  from  0.  These  two  observations  conclude  the  proof. 

Remark  l.The  limit  of  the  above  Cauchy  sequence  is  denoted  by  A(g). 

2. If  G7  is  another  graphical  representation,  then  the  above  theorem  remains  true. 
Moreover,  the  limit  set  does  not  depend  on  the  graphical  representation. 
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3. For  g £ 1ZC(X,Y),  g ^ 0,  one  has 

A(g)=  (J  A(gij)  + {i,j), 

{( hi ) I Si,3#o} 

where  A(gi  j)  + ( i,j ) is  the  translated  limit  set  of  the  graphical  representation  of 

g(X,Y)=gitj. 

Example  l.Let  7?.  = Z2  and  let  the  3-substitution  be  given  by  the  substitution 
polynomial  P^(X)  = id  4-  id X2,  where  id  denotes  the  1 in  Abb(7?.).  The  limit 
set  .4(1),  i.e.  the  limit  of  the  properly  rescaled  graphical  representation  of  the 
sequence  ■ • • 0001000  • • • under  the  substitution  0 000,  1 i -t  101  is  the  triadic 

Cantor  set. 

2. Let  1Z  = Z2  and  define  the  (2  x 2)-substitution  by  the  substitution  polynomial 
Pz(X,Y)  = id  + id  Y + id XY.  The  limit  set  .4(1),  i.e.  the  limit  of  the  properly 
rescaled  graphical  representation  of  the  sequence  g(X,Y)  = 1 (a  single  value  1 at 

(i,j)  = (0,0))  under  the  substitution  0 >-4  q q , 1 ^ ^ , is  the  Sierpinski  triangle. 


3 Cellular  automata 

Substitutions  transform  a sequence  into  another  sequence.  So  do  other  mechanisms 
known  as  cellular  automata.  We  review  briefly  some  essentials  which  we  need 
further. 

Let  S denote  a commutative  ring  with  0 and  1.  According  to  Hedlund10,  a 
cellular  automaton  with  states  in  the  ring  S (actually  it  is  only  necessary  that  S 
is  a finite  set)  is  a continuous  map  A : S( X)  ->•  <S(X)  which  commutes  with  the 
shift  map.  This  implies  that  every  cellular  automaton  A is  given  by  a local  rule 
(j> : SN i+Ni  -»  S,  where  A4 , N2  € N and  A is  defined  as 

m = E fiifj-Ny , fj~Nj+ 1 , ■ • • , fj+Ni-l,  fj+N2)Xj. 

j€  z 

This  means  that  the  j-th  element  of  A(f)  is  obtained  from  / as  a function  q i 
of  the  /-elements  in  a fixed  finite  neighbourhood  of  j.  For  a cellular  automaton 
A : S(X)  — > S(X)  and  an  initial  configuration  / S S(X)  we  define  the  orbit  CLi(/ ) 
of  / (w.r.t.  A ) as  the  two-dimensional  sequence  defined  on  Z x N,  where  row  t £ N 
displays  the  t-th  iterate  of  A on  /.  Formally,  and  in  Laurent  series  notation: 

OO 

0A(l)(X,Y)  = J2At(f(X))Yt- 

t=o 

In  this  note  we  are  dealing  with  cellular  automata  A = Aq,  with  local  rule  <fi  : SD  -» 
S such  that 

A+{f){X)  = '52<Kfi-D+ u...,fj)XP 
j'ez 

Of  particular  interest  are  linear  cellular  automata1,  i.e.,  automata  with  a 5-linear 
local  rule  (j>  ■ SD  S.  In  this  case  the  cellular  automaton  is  also  described  by  the 
multiplication  with  a ’’local-rule” -polynomial  R(X),  i.e.,  A(f,(f)  = R(X)f(X). 
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We  are  particularly  interested  in  cellular  automata  with  the  additional  property 
of  being  k- Fermat9 . They  provide  a framework  to  explain  the  existence  of  fractals 
which  are  generated  by  cellular  automata. 

For  k £ N \ {0}  we  define  the  k-th  power  itk  : S(X)  -4  S(X)  as  7 r/.(f)(X)  = 
£ fj  Xjk  = f(Xk),  i.e.,  7 rk  “inflates”  / by  inserting  (k- 1)  0’s  between  the  elements 
of  the  sequence  / (we  could  call  this  “^-inflation”). 

We  now  recall  the  concept  of  a cellular  automaton  having  the  fc-Fermat  property. 
This  implies  that  row  k of  the  orbit  of  a “^-inflated”  initial  sequence  equals  the 
“^-inflation”  of  row  1 of  the  orbit  of  the  “non-inflated”  initial  sequence,  no  matter 
what  the  initial  sequence  might  be.  Formally: 

Definition  3.1  1,9  A cellular  automaton  A : S(X)  -4  S(X)  is  called  /c-Fermat  if 
(. Ak  o 7 Tfc)(/)  = (71 rjfe  o A){f)  holds  for  all  f £ S{X). 

Remark  1.  If  A is  a /c-Fermat  cellular  automaton,  then  the  sequence  0 is  a fixed 
point  of  A.  Therefore  A(SC(X))  C <SC(X)  for  a fc-Fermat  cellular  automaton. 

2. If  A is  a linear  cellular  automaton  defined  by  the  Laurent  polynomial  R(X), 
then  A is  /c-Fermat  if  and  only  if  R(X)k  = R(Xk). 

Example  l.Let  S = Fp«  , the  field  with  pa  elements  and  of  characteristic  p (p 
being  a prime  number).  Then  any  linear  cellular  automaton  is  p“-Fermat. 

2.  Let  <S  = Zp»,  p a prime  number.  A linear  cellular  automaton  is  p-Fermat  if  its 
local  rule  is  a Laurent  polynomial  R(X)  of  the  form  R(X)  = Q(X)P  , where 
Q(X)  is  a Laurent  polynomial. 

3.  Let  S = Z6,  then  R(X)  = 3 + 4X  defines  a 2-Fermat  linear  cellular  automaton 
and  R(X)  = 2 + 4X  defines  a 3-Fermat  linear  cellular  automaton. 

4.  Let  S = Z5.  Then  the  local  rule  0 : Z§  -4  Z5  defined  by  the  table  with  entries 


x\y 

0 

1 

2 

3 

4 

0 

0 

1 

2 

3 

J 

1 

1 

0 

4 

2 

3 

2 

2 

3 

0 

4 

1 

3 

3 

4 

1 

0 

2 

4 

4 

2 

3 

1 

0 

induces  a (nonlinear)  2-Fermat  cellular  automaton.  In  fact,  it  can  be  shown  that 
there  are  512  different  local  rules  0 : Z§  — > Z5  which  define  a 2-Fermat  cellular 
automaton  on  Z5(X).  Which  demonstrates  that  not  all  /c-fermat  cellular  automata 
have  to  be  linear  cellular  automata. 

We  now  present  a generic  example  which  illustrates  a general  property  of  k- 
Fermat  cellular  automata.  Consider  the  cellular  automaton  with  S = Z3  and  local 
rule  R{X)  = 1 -f  2X  + X2  £ Za[X].  It  is  3-Fermat.  Part  of  the  orbit  Oa(  1)  is 
shown  below: 
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row  00  ■ • • 0001000000000000000000000000000000 • 
row  01  • • ■ 0001210000000000000000000000000000 • 
row  02  • • • 0001101100000000000000000000000000 • 
row  03  • • • 0001002001000000000000000000000000 ■ 
row  04  • • • 0001212121210000000000000000000000 • 
row  05  • • • 0001100000001100000000000000000000 • 
row  06  • • • 0001001000001001000000000000000000 • 
row  07  • • • 0001211210001211210000000000000000 • 
row  08  • • • 0001102201101102201100000000000000 • 
row  09  ■ • • 0001000000002000000001000000000000 • 
row  10  • • • 0001210000002420000001210000000000 • 
row  11  • • • 0001101100002202200001101100000000 • 
row  12  • ■ • 0001002001002001002001002001000000 • 
row  13  • • • 0001212121212121212121212121210000 • 


Row  3 n + l,  l £ {0, 1, 2}  in  this  orbit,  i.e.  the  sequence  T3n+((1),  can  be  obtained 
from  row  3n  as  Al(AMl(\)).  By  the  3-Fermat-property,  row  3n  is  a “3-inflated” 
version  of  row  n (say  ■ • ■ zabcd-  • •),  as  given  by  the  top  sequence  in  the  following 
scheme: 


c 0 0 d 0 0 • • • 

c 2c  c d 2 d d ■ ■ ■ 

b + cb  + c0c  + dc  + d0--- 

(1) 

Also  represented  in  this  scheme  are  rows  3n  + 1 and  3n  + 2,  obtained  by  application 
of  the  rule  R(X).  Observe  from  this  scheme  that  knowledge  of  two  successive 
elements  in  row  n,  say  (a,  b),  is  sufficient  to  determine  all  highlighted  values.  Now 
define  the  substitution  Vl>  which  replaces  the  pair  (a,  b)  by  a (3  x 3)-  array  of  pairs, 
as  follows: 


row  3n 
row  3n  + 1 • ■ 
row  3n  + 2 - ■ 


OOz  0 0 a 0 0 

2 z z a 2a  a 

a + z a + z 0 


b 0 0 

b 2b  b 

a + b a + b 0 


(a,  b)  i-> 


(0,6)  (6,0)  (0,0) 

(a,  6)  (6,26)  (26,6) 

(0,  a + 6)  (a  + 6,  a + 6)  (a  + 6, 0) 


The  last  components  of  the  pairs  in  this  array  are  given  by  the  elements  of  the 
3x3-  array  indicated  in  the  above  three  sequences  (1);  the  first  components  are  the 
elements  that  immediately  precede  them.  It  is  now  clear  that,  if  we  replace  each 
element  (say  6)  in  row  n by  the  corresponding  last  components  in  row  l (l  = 0, 1, 2) 
of  the  substitution  array,  we  produce  row  3 n + l.  In  that  way,  the  orbit  above 
corresponds  to  the  last  elements  of  the  pairs  in  the  arrays  obtained  by  iterating  the 
substitution  starting  from  the  sequence  of  pairs  ■ • • (0, 0)  (0, 1)  (1, 0)  (0, 0)  ■ • • (i.e., 
row  0 rewritten  as  pairs  in  which  the  last  components  correspond  to  the  elements 
of  row  0,  and  the  first  components  to  the  corresponding  preceeding  elements).  We 
show  the  first  step  of  the  iteration  (elements  of  the  actual  orbit  are  underlined) : 


• • • (o,o)(o,  i)(i,o)(o,o)  • • • 
4 
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(0,0)  (0,0)  (0,0) 
(0,0)  (0,0)  (0,0) 
(0,0)  (0,0)  (0,0) 


(0,1)  (1,0)  (0,0) 
(0,1)(1,2)(2,1) 
(0, 1)(1,1)(1,0) 


(0,0)  (0,0)  (0,0) 
(1,0)(0,0)(0,0) 
(0,1)(1,D(1,Q) 


(0,0)  (0,0)  (0,0) 
(0,0)  (0,0)  (0,0) 
(0,0)  (0,0)  (0,0) 


In  the  next  step,  each  pair  is  again  subjected  to  the  substitution  'P.  In  general,  the 
pairs  will  be  D-tuples,  and  the  substitution  will  replace  D-tuples  by  (k  x fc)-arrays 
of  D-tuples  (for  a A;- Fermat  cellular  automaton) . We  now  express  this  property  as 
a formal  theorem. 

Let  D be  a natural  number,  D > 1.  We  define  an  embedding  l = in  of  S(X)  into 
SD(X)  as  follows  (SD(X)  can  be  considered  as  the  set  of  sequences  with  D-tuples 
as  values).  For  f(X)  = € <S(X)  we  define  t(/)  £ SD(X)  as 

t(mx)  = ^(fj.D+1,...,fj)xj. 

je  z 


We  also  define  a projection  p = po  '■  SD(X)  -4  S(X)  such  that  p o i(/)  = / by 
setting  p (Ejez(sh-D+i-  Sj-D+i ,...,sjfi)X^  = T,jeZsjfiXF 

With  these  notions  the  work  of  several  authors7,8,9, 13,14  on  fractals  generated 
by  cellular  automata  can  be  summarized  as 

Theorem  3.2  If  A : S(X)  -4  S(X)  is  a k-Fermat  cellular  automaton,  then  there 
exists  a D £ N and  a regular  ( k x k)- substitution  if  = : SD(X,Y)  -4 

SD(X,  Y)  such  that  for  all  n £ N and  l £ [fc]  we  have 

Ank+l(l)  = p(oMAnm), 


where  oi  is  the  k- substitution  (ipi,i)i£[k]- 

The  interpretation  in  terms  of  fractals  generated  by  cellular  automata  is  clear. 
If  A is  a fc-Fermat  cellular  automaton,  then  the  initial  configuration  1 generates 
an  orbit.  If  the  orbit  is  visualized  by  a graphical  representation,  one  observes  a 
certain  pattern.  According  to  Theorem  3.2,  the  pattern  can  be  described  by  a 
(k  x fc)-substitution  $ = and  by  Theorem  2.4  the  pattern  is  represented  by  a 

compact  subset.  Moreover,  for  linear  fc-Fermat  cellular  automata  it  is  known  that 
the  limit  set  does  not  dependent  on  the  initial  configuration  / 6 iSc(W). 


4 Cellular  automata,  symmetries  and  invariants 

It  is , worthwhile  to  study  cellular  automata  with  a local  rule  <f>  : S2  —>  S which 
exhibits  certain  symmetries  3,4 . 

Definition  4.1  The  map  <f>  : S2  -4  S is  called  rotat.ionally  symmetric  if 
4>((f>(x,y),  x)  = y holds  for  all  x,  y £ S.  If  (f>  is  rotationally  symmetric  and  commu- 
tative, i.e.,  4>(x,y)  — 4>(y,x),  then  <f)  is  called  totally  symmetric. 
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The  following  geometric  interpretation  clarifies  the  definition.  We  represent  the  re- 
lation 2 = y)  as  a triangular  array  xzv,  and  call  it  an  elementary  (^-configuration. 
In  geometric  terms  the  rotational  symmetry  of  (/>  can  be  phrased  as:  Any  rotation 
(by  of  an  elementary  (/(-configuration,  changing  V intoVor  V,  is  an  elementary 
(^-configuration . 

The  commutativity  of  <f>  implies  the  invariance  of  elementary  (/(-configurations 
under  vertical  reflection.  Combined  with  the  rotational  symmetry  of  </>  we  can  state 
that  (j)  is  totally  symmetric  if  and  only  if  the  set  of  elementary  (/(-configurations  is 
invariant  under  the  symmetry  group  of  the  equilateral  triangle. 

Example  l.The  local  rule  </>  : l\  -f  Z2  defined  as  <f>(x, y)  = x + y is  totally 
symmetric. 

2.  Let  Et  = {0, 1,  1 + C}  denote  the  field  with  4 elements.  The  local  rule  i/> : F4  — > 

F4  defined  as  <p(x,  y)  = C,x  + (1  + ()y  is  rotationally  symmetric  but  not  totally 
symmetric. 

3.  There  exist  rotationally  symmetric  local  rules  which  are  not  related  to  linear 
rules.4 

Definition  4.2  A cf>- configuration  of  size  N is  a top  down  equilateral  triangular 
array  with  N elements  in  the  top  row  and  with  values  in  S such  that  any  subtriangle 
xz  is  an  elementary  cf-triangle. 

A </>- configuration  of  size  N is  rotationally  symmetric  if  the  rotated  </>- 
configuration  is  equal  to  the  <j>- configuration . 

A <j>- configuration  of  size  N is  totally  symmetric  if  it  remains  unchanged  under 
the  symmetry  group  of  the  equilateral  triangle. 

Figure  1 shows  examples  of  (/(-configurations. 


Figure  1.  Examples  of  (/(-configurations  of  size  21  for  i />(x,y)  = x + y defined  in  Z2. 

(a):  a nonsymmetric  one,  (b):  a rotationally  symmetric  one,  (c):  a totally  symmetric  one. 
(A  black  cell  represents  the  value  1,  a white  cell  the  value  0). 


Theorem  4.3  4 //(/>:  <S2  — > S is  rotationally  symmetric  and  of  the  form  <j>(x,  y ) = 
rx  + sy,  then  there  exist  for  any  N € N a rotationally  symmetric  <f>- configuration 
of  size  N . 

If  4>  '■  S2  S is  totally  symmetric  and  of  the  form  <p(x,  y)  = rx  + sy,  then  there 
exists  for  any  n E N a totally  symmetric  </>- configuration  of  size  N. 
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As  a next  step  we  introduce  invariants  of  cellular  automata.  The  motivating 
example  for  invariants  is  given  by  the  fc- Fermat  cellular  automata.  Note  that  the 
map  7 r;t  : S(X)  — > S(X)  is  a special  kind  of  ^-substitution.  In  fact,  the  substitution 
polynomial  associated  with  7rk  is  P1Xk  ( X ) = id  and  Ak  oixk  = -j:ko  A.  It  is  therefore 
natural  to  define 

Definition  4.4  Let  A : S — > <S  be  a cellular  automaton.  A k-substitution  £ is  called 
/.■-invariant  of  A if  Ak  o£  = £o  A.  The  set  of  all  k-invariants  of  the  cellular  automaton 
A is  denoted  by  Inv*  (A) . The  substitution  f is  called  a regular  /c-invariant  of  A if  ^ 
is  a k-invariant  and  regular.  The  set  of  regular  k-invariants  is  denoted  by  Inv£(A). 

Remark  If  k = 1,  then  the  1-invariants  of  a cellular  automaton  A with  local 
rule  (f  : SD  ->  S are  given  by  Invi(A)  = {£  : S — > S \ cf>(({ri), . . . , C(ro))  = 
C ru...,rD  £ <S}. 

For  linear  k-Fermat  cellular  automata  the  coefficients  of  the  substitution  poly- 
nomial of  a fc-invariant  £ can  be  characterized  via  the  local  rule. 

Theorem  4.5  If  A is  a linear  k-Fermat  cellular  automaton,  then  Invfc(A)  = 

nLwA). 

Proof.  Let  / £ S(X)  and  let  R(X)  — Li^be  the  local  rule  of  A.  Since  £ 
is  a k-invariant,  we  have  that  R(X)k^(f)  — £(R(X)f(X)).  Since  R is  /c-Fermat, 
the  left  side  becomes  R(Xk)£(f)  = R(Xk)  A’-7(/)(X*)  (using  the  definition 

of  4).  Thus  the  coefficient  of  Xkl+j ,j  £ [fc],  in  R(Xk)£(f)  can  be  computed  as 
rotj(fi)  + + • • • + rD£j(fi~D).  On  the  other  hand,  the  same  coefficient 

in  {(R(X)f(X))  turns  out  to  be  £j(r0/;  + • • • + roft-D),  he.  & o A = A o 
Therefore,  £ is  a /c-invariant  if  and  only  if  each  £,•  is  a 1-invariant  of  A.  • 

Example  l.For  S — TLi  and  R(X)  = 1 + X,  we  have  Imq(A)  = {s  h4  s,  m 0} 

( s £ S ) and  Invi(A)  = InvJ(A).  The  automaton  specified  by  this  rule  is  2"-Fermat 
(for  any  n £ N).  Thus  Inv!j(A)  contains  22  regular  2-invariant  substitutions  : 
0 M-  00,  1h>  00;  0 i t 00,  1 10;  0 •->  00,  01  and  0 >->  00,  1m-  11.  InvJ(A) 

contains  24  regular  4-invariants,  among  which  0 i->  0000,  1 1011. 

2.  For  <S  = Z3  and  R{X)  = 2 + 2X,  we  have  Invj(A)  = {s  as  + 0 \ a,  0 £ Z3} 

and  Inv^(A)  = {s  h*  as  \ a £ Z3}.  Inv3(A)  contains  33  regular  /c-invariants, 
among  which  0 000,  1 i->  201,  2 i-»  102. 

3 . R(X)  = 1 + 3X  £ Ze[X\  defines  a 2-Fermat  automaton  such  that  |Invi(A)|  = 54 
and  |Invj(A)|  = 18.  Thus  there  are  218  regular  2-invariants. 

In  analogy  with  the  substitution  T introduced  before  Theorem  3.2,  we  associate 
here  with  a linear  fc-Fermat  automaton  A,  given  by  R(X)  — s + rX,  and  a k- 
invariant  £ £ Inv£(A),  a (lx  /c)-substitution  = 4>(A,£)  : S2(X,  Y)  — ► S2(X,  Y). 
Let  (a,  b)  £ S 2 and  consider  the  polynomials  Qm(X)  = R(X)m£(aX~1  + b)  = 
Yjjgz  1j,mX3  for  m £ [k],  For  m £ [fc]  the  coefficients  q-i,m,q0,m, . . . ,qk-i,m  of 
Qm(X ) are  uniquely  determined  by  a and  b , thus  these  coefficients  can  be  regarded 
as  functions  from  S 2 to  S.  Therefore,  the  substitution  = (tAm)(,mG[fc].  defined 
as 


V’i.m  ; >52  -*  S2 

(a,b)  h->  (qi-i,m(a,b),qitm(a,b)) 
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for  l,m  £ [&],  is  well  defined  and  called  the  induced  substitution. 


Remark  1.  If  £ is  a regular  fc-invariant,  then  the  induced  substitution  '5^  is  a 
regular  substitution. 

2.  If  the  local  rule  of  a linear  cellular  automaton  is  given  by  the  polynomial  R(X)  6 
5[X],  of  degree  d,  and  if  £ is  a regular  ^-invariant,  then  the  induced  substitution 
is  defined  by  maps  from  <Sd+1  which  are  constructed  in  a similar  manner  as  for  the 
case  d — 1 considered  above,  (d  + l)-tuples  are  then  replaced  by  a (k  x fc)-arrav 
of  (d  + l)-tuples. 

3. If  A ^ is  any  cellular  automaton  such  that 


and  with  a ^-invariant  £,  then  the  above  construction  yields  an  induced  substitu- 
tion as  well. 

Example  For  S = Z2,  the  polynomial  R(X)  = 1 + X defines  a 22-Fermat  linear 
cellular  automaton  A.  By  Example  1 following  Theorem  4.5,  we  have  that  the 
substitution  £ defined  by  the  polynomial  P^(X)  = id  + OX  + idX2  +idX3  defines  a 
4-invariant  of  A.  Consider  a sequence  • • • zabcd  ■ ■ •,  then  substitution  £ transforms 
this  into  the  top  row  shown  below: 


•••  z 0 z z a On  a 
zzOz+aaaO 
Ozz+azOa 
z a a z a 


b 0 b b 
ja  + bbbO 
a -I-  b a 0 b 
b b a b 


b 

b 


c 0 c c • • • 
+ ccc0'-- 
+ c b 0 c ■ ■ ■ 
c c b c ■ ■ ■ 


The  other  rows  are  part  of  the  orbit  generated  with  the  given  cellular  automaton 
rule  starting  from  the  first  row.  The  induced  substitution  '1^  can  be  read  directly 
from  these  rows: 


(a,  b)  i-> 


(a,  b)  (6,0)  (0,6)  (6,6) 

(0,  a + 6)  (0  + 6, 6)  (6, 6)  (6, 0) 
(a,  a + 6)  (a  + 6,  a)  (a,  0)  (0, 6) 
(a,  6)  (6,6)  (6,  a)  (a,  6) 


(2) 


The  following  theorem  is  a straightforward  generalization  of  Theorem  3.2.  It 
exploits  the  existence  of  ^-invariants  instead  of  the  fc-Fermat  property. 

Theorem  4.6  Let  A#  be  a cellular  automaton  and  let£  be  a regular  k -invariant  and 
let  'Ilf;  = (if>i,j)ilje[k]  be  the  induced  k-substitution.  There  exists  a natural  number 
D,  and  embedding  l = ld,  such  that  for  / £ S(X)  with  £(/)  = / and  all  n £ N, 
l £ [£]  we  have 


Ank+l(l)=p{aMAn(m), 

where  07  is  the  k-substitution  (V’i.i )ie[*]  • 

Proof.  Since  £ is  a fc-invariant  and  / = £(/),  we  have  Akn(f ) = Akn(£(f))  = 
£(An(f))  for  all  n £ N.  That  proves  the  assertion  for  l = 0.  The  assertion  for  l / 0 
follows  from  the  definition  of  the  induced  substitution.  • 

We  illustate  this  theorem  by  continuing  the  above  example,  with  <S  = Z2, 
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R(X)  = 1 + X and  the  4-invariant  which  produces  the  following  substitutions: 

0 H*  0000,  1 1-4  1011.  Take  / as  the  sequence  which  is  the  limit  obtained  from 
applying  this  substitution  to  • • • 0001011000  • • • (the  boldface  element  is  at  position 
0).  This  yields  the  top  row  in 

/ • ■ • 0001011000010111011000000000000000010111011 • • • 

A(f)  ■ ■ ■ 0001110100011100110100000000000000011100110 • • • 

Observe  that  f = £(/).  Theorem  4.6  states  that,  for  example,  row  6 (=  1 • 4 + 2) 
in  the  orbit  OaU)  can  be  obtained  from  A(f_)  displayed  above  , by  imbedding  this 
last  sequence  in  a sequence  of  pairs  (corresponding  to  D = 2),  as  given  by 

• • • (0, 1)(1, 1)(1, 1)(1,Q)(0, 1)(1,0)(0,0)(0,0)(0, 1)(1, 1)(1, 1) 

(the  last  (underlined)  components  in  these  pairs  are  the  elements  of  A(f),  the 
first  components  are  the  ones  that  immediately  precede  them).  Then  apply  the 
substitutions  displayed  in  row  2 of  the  induced  substitution  array  4>?  given  above 
(see(2):  row  2 is  actually  the  third  row,  as  counting  starts  at  zero).  This  gives 

•••(0,1)(1,Q)(0,0)(0, 1)|(l, 0)(0,1)(1,0)(1, 1)|(1, 0)(0,I)(1,0)(1, 1)|(l, i)(l, I)  ••• 

The  last  components  in  these  pairs  form  row  6 in  Oa{[)-  In  a similar  way,  using  row 

1 of  the  ^-substitution  array,  where  l E {0, 1,2,3),  starting  from  the  embedding 
of  row  An(f)  in  a sequence  of  pairs,  would  produce  row  4 • n + l in  this  orbit. 

As  a consequence  of  Theorem  2.4,  we  obtain 

Corollary  4.7  //£  is  a regular  k-invariant  of  the  cellular  automaton  A and  if 
is  the  induced  substitution,  then  the  sequence  ^prG(4>£  (t(/)))^  ^ is  a Cauchy 
sequence  for  all  / E SC(X). 

As  a final  step,  we  combine  the  notion  of  symmetry  of  a cellular  automaton  with 
the  existence  of  induced  substitutions.  We  are  interested  in  the  following  problem. 
Suppose  /o, . . . , /iv-i  is  the  top  row  of  a (^-configuration  of  size  N and  assume  that 
the  configuration  is,  e.g.,  rotationally  symmetric.  If  £ is  a regular  ^-invariant  of 
A$,  then  £()T)  fjXi)  gives  anther  top  row  of  a ^configuration  of  size  kN;  what  are 
the  symmetry  properties  of  the  larger  configuration? 

Definition  4.8  Let  A$  be  a cellular  automaton  with  local  rule  (j> : S2  -4  S and  let 
the  substitution  £ = (£;);e[fc]  be  a k-invariant  of  A.  The  k-invariant  is  palindromic 
if  £;  = £fc_j_ i holds  for  all  i E [fc]. 

The  k-invariant  £ is  called  rotationally  symmetric  if  the  induced  substitution 
satisfies 

p(V»i,o(a,  b))  = pW'/t-i  ,i(a,b)) 
for  all  (a,  b)  E S 2 and  all  i E [fc]. 

If  £ E Lav*,  (A)  is  palindromic  and  rotationally  symmetric,  then  £ is  called  totally 
symmetric. 

An  example  of  a rotationally  symmetric  4-invariant  is  given  by  the  example  pre- 
ceeding  Theorem  4.6:  the  last  components  of  the  rightmost  column  in  the  4^ -array 
equals  the  last  components  of  its  top  row. 


62 


Remark  If  £ = (£;)  is  a ^-invariant  for  A$  and  = (ijji,m)  the  induced  sub- 
stitution, then  we  have  p(ipit0(a,  b))  = fi(b)  for  i £ [fc]  and  there  exist  functions 
Pi  £ Abb(<S),  i £ [fc],  such  that  p(^_lii(a, b))  = Pi(b).  In  particular,  if  £ is  rota- 
tionally  symmetric,  then  we  have  pi  = £*  for  all  i £ [k]. 

Lemma  4.9  4 Let  A be  a cellular  automaton  with  local  rule  f : <S2  — > S.  If 
f(X ) G <S[X]  is  of  degree  < N,  then  f(X)  defines  a rotationally  symmetric  f- 
configuration  of  size  (N  + 1)  if  and  only  if  for  the  coefficients  Oij  of  0A(f)(X,  Y ) 
the  equations  Oifi  = Ona  hold  for  i = 0, . . . ,N . 

In  the  linear  case  a carefully  performed  substitution  can  preserve  symmetries. 

Lemma  4.10  Let  be  a cellular  automaton  with  rotationally  symmetric  local 
rule  f(x,y)  = rx  + sy.  If  fo,  ■ ■ ■ , fN-i  is  a top  row  of  a rotationally  symmet- 
ric (}>- configuration  of  size  N and  if  £ £ InvJ/A^)  is  rotationally  symmetric,  then 

9o,---,9kN-i,  where  g{X)  = J2j<=z9jXj  = £(£^0 1 fiX1),  is  the  top  row  of  a 
rotationally  symmetric  <j>- configuration  of  size  kN. 

Proof.  Let  f(X)  = XjL"1  fix*  and  O = 0A  {f)(X,Y)  = £ OuXlY^  its  orbit.  If 
denotes  the  induced  ( k x /^-substitution,  then  'I'^(O)  = £ QijXjY^  is  the  orbit 
of  £(/).  It  remains  to  prove  that  0iiO  = QkN-i, i holds  for  all  i = 0, . . . , kN  - 1. 

By  Lemma  4.9,  we  have  Oifi  — On- i,i  for  i = 0, . . . , N - 1.  Since  £ is  rota- 
tionally symmetric,  we  have  &kj+i,o  = £f(Oj,o)  for  kj  + l £ [Nk]  and  l £ [k]  and 
QkN-i,kj+i  — ^i(On-ij)  = fi(Oj:o),  which  proves  the  assertion.  • 

Lemma  4.10  remains  true  if  rotationally  symmetric  is  replaced  by  totally  symmetric. 

We  can  now  establish  the  existence  of  fractals  with  prescribed  symmetries. 

Theorem  4.11  Let  A with  <f>(x,y)  — rx  + sy,  define  a cellular  automaton  with 
a rotationally  (totally)  symmetric  linear  local  rule.  Let  f(X)  £ S[X]  of  degree 
< N define  a rotationally  (totally)  symmetric  cj>- configuration  of  size  (N  + 1) 
and  let  £ £ InvJ/A^)  be  a regular  rotationally  (totally)  symmetric  k-invariant. 
If  i(0A(f))  = i>2  (0A(f))  denotes  the  embedding  of  the  orbit  0A(f)  which  replaces 
( 0A(f))j,i  by  (( 0A(f))j-i,i , ( 0A(f))j}i ),  and  if^>^  denotes  the  induced  substitution, 
then  the  sequence 

(fWUi  G^g(*<”m<u(/)>>>' n v)n6N, , 

where  V is  the  equilateral  triangle  given  by  the  points  (0,0),  (1,0)  and  (1/2,  a/3/2) 
(ordinate  axis  positively  oriented  downwards),  is  a Cauchy  sequence.  The  limit, 
denoted  as  X = A(/,£),  is  a compact  set  in  V and  X is  rotationally  (totally) 
symmetric. 

Proof.  By  Theorem  2.4,  the  limit  1 /knG(^(i(f)))  exists  and  is  contained  in 
[0,  N + l]2.  The  factor  transforms  the  triangle  {(a :,y)  £ [0,N  + l]2  \ y < x}  into 
an  equilateral  triangle  of  size  1.  Since  the  limit  is  independent  of  the  particular 
graphical  representation  (“pixel-shape”),  we  can  define  7 : S2  \ {0}  ->•  TL(¥f)  by 
setting  7(s_i,s0)  = U{(/0)|sj  / 0;j  = -1,0}.  Then  the  rescaled  version  of 
G7(*cnW/)))  n V is  rotationally  (totally)  symmetric  for  each  n £ N.  Thus  the 
limit  is  rotationally  (totally)  symmetric.  • 
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Figure  2.  Symmetric  fractals  generated  by  cellular  automata,  (a)  and  (b):  states  in 
Z2  = {0,1}  with  local  rule  (f>(x,  y)  — x + y and,  for  (a):  with  top  row  gener- 

ated from  the  initial  top-row  configuration  1101101101  and  the  16-substitution  0 
0000000000000000,  1 t-4  0110011000000110,:  for  (b):  initial  configuration  0110110110 
and  8-substitution  0 h->  00000000,  1 i-4  01111110.  (c):  states  in  F4  = {0, 1,C,  1 + C}  = 
{0, 1,  2,3}  with  local  rule  <f)(x,y ) = (x  + (1  + Qy,  initial  state  0130  and  4-substitution 
0 0000,  1 i-t  1021,  2 f y 2032,  3 t-J-  3013.  (d):  states  in  Z3  = {0,1,2} 

with  local  rule  <f>(x,y)  = 2x  + 2 y,  initial  configuration  11011011  and  9-substitution 
0 h-»  000000000,  1 002020200,  2 i->  001010100.  (a)  and  (c)  are  rotationally  sym- 

metric, (b)  and  (d)  are  totally  symmetric.  The  Hausdorff  dimensions  are  respectively 
log  117/  log  16  « 1.717,  log  36/  log  8 « 1.723,  log  13/  log  4 ss  1.850,  log  54/  log  9 « 1.815. 


Figure  2 shows  a few  examples  of  symmetric  fractal  limits  obtained  from  rescaled 
cellular  automata  orbits  generated  by  a proper  initial  top-row  configuration  sub- 
jected to  a proper  substitution  as  presented  in  the  above  theorem. 
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5 Conclusions 

The  above  described  method  allows  one  to  construct  fractals  which  reflect  the  local 
rule  of  a cellular  automaton  as  well  as  showing  a global  symmetry.  Since  these 
fractals  are  generated  by  a (k  x ^-substitution,  these  fractals  can  also  be  described 
in  terms  of  hierarchical  iterated  function  systems,  in  particular  by  so  called  fc-adic 
HIFS8.  By  introducing  a transition  matrix8  for  the  generating  substitution  it  is 
possible  to  compute  the  fractal’s  box-counting  dimension5’6  which  coincides  with 
the  Hausdorff  dimension. 
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